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Abstract 

The design of the LHC and future colliders aims their operation with high intensity 
beams, with bunch population, N p , of the order of 10 11 . This is dictated by a 
desire to study very rare processes with maximum data sample. HEP colliders are 
engineering structures of many kilometers in length, whose transverse compactness 
is achieved by the application of the superconducting technologies and limitations 
of cost. However the compactness of the structural elements conceals and potential 
danger for the stable work of the accelerator. This is because a high intensity beam 
of positively charged particles (protons, positrons, ions) creates around itself an 
electric self-field of very high intensity, 10 5 - 10 6 V/ cm. Being located near the 
conducting surfaces, at the distances of 1-20 mm away from them, the field of 
such bunches activates the field emission of electrons from the surface. These 
electrons, in addition to electrons from the ionization of residual gases, secondary 
electrons and electrons knocked out by synchrotron radiation, contribute to the 
development of a dense electron cloud in the transport line. The particles of 
the bunch, being scattered on the dense electron cloud with N e ~ N p , leaves 
the beam and may cause noticeable damage. The paper presents an analysis 
of the conditions, under which the field emission in the LHC collimator system 
may become a serious problem. The analogous analysis of a prototype of the 
International Linear Collider (ILC) project, USLC, reveals that a noticeable field 
emission will accompany positron bunches on their entire path during acceleration. 
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1 Introduction 



Searches for rare processes and the production of new particles require collider experiments 
to be runs with the highest possible luminosity, given by the standard expression 

r n b -N x -N 2 

L ~ . (1) 

a x o y 

To achieve this goal it is necessary to have beams with maximal bunch population, Ni, large 
number of bunches, n&, and bunches with minimal transverse sizes, a x , a y . The electric and 
magnet self-fields created by an intense particle bunch, with modification by the accelerator 
environment (beam pipe, accelerator gaps, magnets, collimators, etc.) due to induced surface 
charges or currents, can have very high values at a distance of the order of a few millimeters 
from the bunch. 

It is a well known fact that under the influence of an electric field F of the order of 
10 5 V/cm at room temperature the smooth surface of solids (conductors, semiconductors, 
dielectrics) and liquids emits electrons 3 (see ^ and references therein). The phenomenon, 
known now as electron field emission (FE), was discovered by Robert Wood j^j in 1897. Later, 
in 1928, Fowler and Nordheim 3 proposed a quantum theory of field emission from metals in 
terms of electrons tunneling through a potential barrier. Application of a high electric field 
to the metal produces a triangular shaped potential energy barrier through which electrons, 
arriving at the metal surface, may quantum mechanically tunnel. 

Such electron emission usually precedes electrical breakdown in vacuum and is called 
pre-breakdown emission. We note that the maximum field strength that can be maintained 
between two conductors in air is limited to less than about 10 4 V/cm, above which dielectric 
breakthrough leads to the formation of a plasma. With semiconductors, of the order of 10 6 
V/cm can be maintained. Fields of the order of 10 8 V/cm = i v/A can also be established 
within 10 3 A of a metal tip with a tip radius of less than 10 3 A, provided breakthrough is 
avoided by working in ultrahigh vacuum. In fields larger than 1 V/A a variety of processes 
start to develop: field desorption 4 , field evaporation 5 , vacuum discharge. 

In this paper, we discuss the field emission phenomena activated in strong fields created by 
a positively charged particle beam. Primary attention is paid to the field emission processes 
and beam induced multipacting in the LHC collimator system. The primary and secondary 
collimators in the IR3, IR6 and IR7 insertions are made of graphite [77j-[S3] with parallel jaw 
surfaces and small gaps, ~2 mm, at 7 TeV. As shown below, for the LHC bunch parameters 
and structural features of the LHC collimator system, the density of the electron emission 
current can reach a significant magnitude, which in turn can lead to a "thrombosis" of the 
collimators. 

The paper is organized as follows. In Section 2 we briefly review the classical Fowler- 
Nordheim theory for point cathodes, following which we give an estimation of the tunneling 
time and discuss deviations from the F-N theory due to extensions to non-zero temperatures 
and "real" broad-area cathodes. In Section 3 we discuss electrical self-fields generated by 
charged particles encountered in high energy physics: an elementary particle, a relativistic 
dipole, a particle bunch. In Section 4, the results of studies made with the use of the modified 
F-N theory are presented. Here we discuss the field emission in the LHC collimator system 



3 We denote the electric field by F follow the established tradition in publications devoted to field emission. 

4 Field desorption is a process in which an atom or molecule leaves a solid surface due to the influence of a high 
electric field. The atom or molecule usually departs from the surface in a charged state as an ion. 

Field evaporation is the removal of lattice atoms as singly or multiply charged ions from a metal in a strong 
electric field F of the order of several V/A, as occurs at field ion tips. 
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and the main linac beam pipe of the US linear collider. Finally, Section 5 discusses the 
dynamics of emitted electrons, the electron multiplication chain, and the development of an 
electron cloud in the self-fields of a bunch. 



2 Field Emission in Strong Electric Fields 
2.1 Fowler-Nordheim theory 

In the framework of Fowler-Nordheim (F-N) theory, the current density of field emission of 
electrons from a metal can be written in the following form |3].[1]-[TT] 

Jfn = ej n(£ x )D(£ x ,F)d£ x , (2) 

where D(£ X ,F) is the penetration coefficient and n(£ x ) is the number of electrons at the 
energy £ x incident in the x-direction on the surface barrier from inside of the metal. 

An electron outside a metal is attracted to the metal as a result of the charge it induces on 
the surface (image force). In the externally applied accelerating electric field F, the potential 
energy of the electron is 6 

e 2 

V(x) = — — — eFx, when x > 0, (3) 
Ax 

where x denotes the distance from the surface and the first term accounts for the image 
potential. With use of the potential energy © and the Fermi energy distribution of electrons 
in the conduction band, one finds [H],[I]-J5 that 

F 2 r V 3/2 ^ 

j fn (f) = ^ ex p{-v ( 4 (4) 

where J is the current density in A/cm 2 , F is electric field on the surface in V/cm, and ip is 
the work function in eV. The field-independent constants A and B and the variable y are 



A = ^—r = 1.5414- 10~ 6 , B = °" v f = 6.8309- 10 7 , y = = 3.7947- 10~ 4 — (5) 

8irh 3en ip ip 

where — e is the charge on the electron, m is the electron mass and h is Planck's constant. 
The numerical values of A and B correspond to recent values of the physical constants |13j . 
We note that under field emission conditions, 0< y <1. 

The Nordheim function takes into account a lowering of the potential barrier due 
to the image potential (the Schottky effect) and its distinction from an idealized triangular 
shape. The function t(y) in the denominator of equation Q is defined as 

t(y) = <%) - (2y/S)(d#/dy). (6) 

The function #(y) varies from i?(0) = 1 to i9(l) = with the increase in field strength, 
however t(y) is quite close to unity at all values of y. 

For a typical metallic <p of 4.5 eV, fields of the order of 10 7 V/cm are needed to have 
measurable emission currents. In considering magnitudes, one must always keep in mind the 
rapid variation of the exponential function. For instance, an increase in F of only a factor of 
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3 In the present paper we adopt the Gaussian CGS system. 
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two from 1 x 10 to 2 x 10 7 V/cm increases the current density by 15 orders of magnitude 
(from 1(T 18 to l(T 3 A/cm 2 ) ! 

At a field strength of the order of F cr = ip 2 /e s = 6.945 • 10 6 • ip 2 V/cm the height of 
the potential barrier vanishes and i?(l) = 0. For instance, for copper ipcu = 4.65 eV giving 
F cr (Cu) = 1.5 • 10 8 V/cm, and similarly for graphite, ip gr = 4.6 eV, F cr (gr) = 1.47 • 10 8 
V/cm. At this field level one would expect the orderly bound states characteristic of the 
solid to lose their integrity. 

For a long time only tabulated values of #(y) and t(y) ^1] were used in calculations, 
see |5]-|11|. Recently |15| . several parameterizations of the functions $(y) and t(y) were 
proposed. 

The theory of field emission from metals has been subjected to fairly extensive verification. 
A variety of methods have been employed over many years for the experimental measurements 
of the emission current as a function of the field strength, the work function and the energy 
distribution of the emitted electrons [7]-|llj. The F-N theory (jlj) of electron emission from 
plane and uniform metal surfaces (single-crystal plane) at T ~ may therefore be considered 
well established on experimental basis as well as on theoretical grounds. 

2.2 Tunneling Time 

Using the Heisenberg uncertainty principle one can estimate the tunneling time. For electrons 
near the Fermi level, the uncertainty in their momenta, Ap, should correspond a barrier of 
height ip where 

ip = l^p Ap = (7) 
If the corresponding uncertainty in position, 

Ax ~ — j (8) 

2^/2^ v ; 



is of the order the barrier width 



A*~^, (9) 



one expects to find electrons being emitted. 

On the other hand, the uncertainty in energy is of the order 

Ae « eFAx (10) 

and therefore the estimated value of the tunneling time is 



A^^p«3.37.10- 8 #. (11) 
eF F 

Here the field strength F is in V/cm, the work function ip is in eV and the tunneling time 
At is given in seconds. As an example, for graphite and a field of 10 6 V/cm this gives 
At = 7.23 • 10 -14 s, quite a short time scale, when the field is large. If the field is created by 
a relativistic bunch of length L, the field is acting on a certain area of the surface during the 
time t = L/c. For a LHC bunch of L=7.55 cm one obtains r = 2.52 • 10~ 10 s, which is long 
in comparison with the emission time At. 



4 



2.3 Deviations from the Fowler-Nordheim theory 



2.3.1 Temperature dependence 

The main equation Q of the F-N theory was derived for an idealized metal in the framework 
of the Sommerfeld model, with an ideally flat surface and at a very low temperature, T ~ 0. 
The temperature dependence of the field emission current (FEC) is completely connected 
with the change of the spectrum of electrons in the metal with an increase in T. Therefore, 
at non-zero temperatures the F-N theory must be modified to take into account the thermal 
excitation of electrons above the Fermi level. For the so-called extended field emission region, 
Murphy and Good [H] (see also obtained the following elegant equation 

J FN (F,T) = -^-J FN (F,0), (12) 

which account for the temperature dependence of the FEC. Here u> = k B T/k B T and 

2 F 
3Bt{y)W<p 



where k B is the Boltzmann's constant, and T is the absolute temperature in K. It can be 
shown 5 j that equation (|12j) is a valid approximation when the following two conditions are 
satisfied: 

l-l-i _ ¥>(l-y) 2/2\i/2 



u < 
and 



1 + ri 



u < 



1 

1 2 



-i / 2cp \ 1/2 

■ r >^vk)) ■ (15) 

At very low temperatures, when ttlo is small, equation (|T2*|) reduces to @. By expanding 
sin7ro; in a series, one gets for practical use the formula 

J FN (F,T)/J FN (F,0) = 1 + 1.40 • 10 8 (^/F 2 )T 2 . (16) 

It is easy to estimate using (fTK|) that for (p =4.5 eV at room temperature, T = 300K, and 
F = 1 x 10 7 V/cm and 2 x 10 7 V/cm, the temperature factor in (J12[) equals 1.57 and 1.14, 
respectively. Thus, the temperature factor appear to be a sizeable correction. 

The energy distribution of the emitted electrons is determined by two effects: the low- 
energy slope by the tunneling probability, and the high-energy slope by the electron distri- 
bution in the metal, and thus by the temperature T of the emitter. For increasing fields, the 
width of the energy distribution will grow approximately with the field, while the position 
of the spectrum stays close to the Fermi level, E F |llj : 



dJ _ J FN (F,0) ex.p(e/k B T ) 



de ek B T Q 1 + exp(e//c fl T) 
Here e = E — E F and E is the kinetic energy of the electron. 



(17) 



2.3.2 Anomalous High FEC 

Since the first experimental verifications of the F-N field emission theory, it was noticed 
that for broad-area cathodes the electrical breakdown typically occurs for field values one 
or two orders of magnitude smaller and pre-breakdown FEC is many orders of magnitude 
greater than the values predicted by the F-N equation. Nevertheless the observed emission 
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current follows the F-N law, provided one makes the substitution F — ► F,, = (3 FN F for all 
occurrences of F in (|12[) """J1I201 1211 122| . Here /? FJV is known as the field enhancement factor. 
Generally, f3 FN values in the range 50 < (3 FN < 1000 have been observed in many experiments. 
Emission does not occur homogeneously over the surface but is rather concentrated in /im and 
sub-^m sized spots [H I20| "24"]. Then electrons are concentrated in tiny jets and with further 
increase in voltage, this ultimately leads to breakdown. Physics of the local field enhancement 
and a vacuum arc discharge, in general, is very complex and still not completely understood. 
For reviews see |20|.|25|. |28j . As one sees from ("""J , FEC increases both with an increase of 
F as well as with decrease of (p, or due to a combination of these two factors. There are a 
number of mechanisms which lead to an increase of FEC. Here we enumerate some of them. 

Geometrical field enhancement. The real surface of solids is not perfect. For many 
years, the anomalous emission was universally explained by invoking a 'projection model' 
["""Jl "2*01 [2*2*1 "2*""] . This model assumes the presence on the cathode surface of a number of 
microscopic projections (crystalline defects, impurities, whiskers or dust particles), sharp 
enough to cause a geometric enhancement of the local field at the projection tip to a value 
some 100 times greater than the nominally applied field. The projection model was confirmed 
by experimental evidence [*Q 1201 "2*""] obtained with shadow/scanning electron microscopy. 
Even on optically polished cathodes made of stainless steel, tungsten, copper or aluminum, 
needle-like projections about 2 /xm high have been found. These projections are capable of 
producing field enhancements of the order of 100 at pre-breakdown emission sites. 

Resonant tunneling. To account for enhanced emission from gas condensation, another 
model assumes that ad-atoms are responsible for creating localized energy levels near the 
surface. One dimensional calculations [S2] show that the tunneling process of electrons with 
energies close to the localized states can be resonantly enhanced. The calculations predict 
that tunneling is enhanced by up to a factor of 10 for adsorbates less than a single monolayer 
thick. For example, adsorbed water with its strong dipole moment can be a critical factor 
in enhanced electron field emission ["""""]. Water is certainly one of the main adsorbates on 
the conductor surface, especially if the surface is not baked. There is also observation that 
adsorbed oxygen enhances field emission |34j . 

Surface states (SS). By the level of complexity, this is a completely different mechanism 
resulting the FEC increase due to a lower electron work function of SS. A semi-infinite crystal 
extending from x = — oo to x = can be viewed as a stack of atomic layers parallel to the 
given crystallographic plane but, in contrast to the infinite crystal, all layers cannot in general 
be identical. Consequently, the potential the electron sees at the metal-vacuum interface is 
unavoidably different from that in the bulk material because the electronic charge distribution 
is different on the surface. On metal surfaces, the density of electrons is high, of the order 
of 10 15 cm -2 . That gives rise to creation and decay of a variety of electronically excited SS 
"SHI HOI in] • The term applies because the wave function of such a state is localized near 
the surface, decaying "exponentially" on either side of it. The surface states are classified 
as Tamm ("""""] and Shockley [HZ| surface states (TSSS), and image potential states (IPS) 
|38| . Tamm states are more localized and arise when the potential in the top surface layer 
is significantly distorted, whereas Shockley states show a strong free-electron-like dispersion 
and may be interpreted as dangling bonds of surface atoms. Such SS are described by 
wave functions whose center of gravity lies in the immediate neighborhood of the metal- 
vacuum interface on the metal side of it. Self-consistent calculations and angle-resolved 
photoemission measurements |42| . have shown that such states exist on the (100) plane of 
copper. Employment of scanning tunneling microscopy at low temperatures allows one to 
detect the surface topology and view SS of many noble metals. For a state-of-the-art review 
see 02]. 
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Image potential states j3H] (see also [lEl-Jin], [SH]), arise for an additional electron in front 
of the surface. The screening of the charge by the metal electrons give rise to a long-range 
Coulomb-like attractive image potential which leads to bound states which form a Rydberg 
series with energies E n , where 



E n = 7 '-^>eV, n = l,2,..., < a < 0.5 (18) 

(n + a) z 

converging toward the vacuum energy. Multiphoton photoemission spectroscopy allows us 
to identify IPS on graphite j5U]. For IPS on graphite a quantum defect of a = —0.04 ± 0.05 
is obtained. This demonstrates that IPS on graphite exhibit an almost ideal hydrogen-like 
behaviour. The maximum of the IPS wave function density with n=l is located several 
Angstroms away from the surface at the vacuum side |471 164j . 

While the energy levels of TSSS are spread around the Fermi level, the energy levels of 
IPS are allocated very close to the vacuum level, E n > —1.0 eV. That is, the effective work 
function of TSSS or IPS is equal to (p e ff = \E n \. 

Fullerenes, nanoparticles, nanotubes etc [Hll- [S3J - The LHC collimator jaws are 
made of graphite |78]-|H4j. Graphite can be polished only mechanically, and therefore may 
serve as a large area emission surface. The simplest explanation for this is that the surface 
morphology [MJ consists of countless fractures, protrusions, edges, and other features which 
are readily demonstrated to emit electrons with field enhancement. Fig ^ shows the typical 
graphite surface at high magnification |H2]- In conjunction with the previous discussion, we 
shall consider briefly some properties of graphite and its derivatives. 

Graphite consists of layers which are formed by regular hexagons. Partial thermal decom- 
position of graphite layers (as a result of heating at higher local emission current densities) 
can produce on the surface fullerene molecules, nanoparticles and also long tiny tubes, nan- 
otubes. Under the influence of a high voltage field these objects are able to emit electrons. 

The Cqq molecule occupies the central place among fullerenes. This molecule resembles 
the commonly-used surface of a football. The carbon atoms are distributed on a spherical 
surface at the vertices of twenty regular hexagons and twelve regular pentagons. The radius 
of the Cqo molecule, deduced from X-ray structure analysis [HE], is 0.357 nm. The common 
structural components found in graphite and the Cqq fullerene molecule determine the nature 
of the process of formation of fullerenes by the decomposition of graphite. Moderate heating 
of graphite breaks the bonds between the separate layers and an evaporated layer then splits 
into separate fragments. These fragments are combinations of hexagons one of which becams 
used to form a cluster. 

The process of formation of fullerenes from graphite yields further various structures 
which are composed, like graphite, of six-member carbon rings. These structures are closed 
and empty inside. They include nanoparticles and nanotubes. The nanoparticles are closed 
structures similar to fullerenes, but of much larger size. In contrast to fullerenes, they may 
contain several layers. Such multilayer spherical structures are called onions. The study 
|55j of field emission properties of diamond-like carbon films indicate that the films have a 
faceted morphology with carbon cluster size in the range 150-400 nm. The fitting of I-F 
(current-electric field) data to the F-N equation gives j3 FN ~160, assuming if = 5 eV. 

Carbon nanotubes (CNT) are elongated structures with surfaces formed by regular hexagons. 
Separate extended graphite fragments, which are then vent into a tube, form the base of CNT. 
These tubes are up to several microns long with diameters of a few nanometers. They are 
multilayer structures with rounded ends. One end is attached to a surface and the other is 
free. 
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Figure 1: Typical rough graphite surface |62j . 



As follows from direct measurements |57|-|63|. CNT have quite a high field emission 
capability. This property centers around the extraordinarily small diameter of CNT, which 
provides a sizeable enhancement of the electric field strength near the CNT cap in relation to 
that averaged over the entire volume of the interelectrode gap. For instance, at a voltage near 
650 V, an emitting area of about 1 mm 2 provides an emission current of order 0.5 mA. The 
measured emission characteristics correlate well with the F-N equation, with the condition 
that the value of the electric field strength is taken at the point of the electron emission. 
Since this point is situated close to the sharpened top of the CNT cap, the local magnitude 
of the electric field considerably exceeds its mean value. Thus, the above mentioned effect 
of field enhancement becomes involved, with the magnitude of (5 FN estimated to be of the 
order of 1000, assuming the work function to be close to <p = 5 eV. 

High emission properties of CNT cannot be attributed only to their high aspect ratio. 
As demonstrated in model calculations and measurements (641 165| , single- and multiwalled 
CNT generate both IPS as discussed above and a new class of surface states, the tubular 
image states, owing to their quantized centrifugal motion. Measurements of binding energies 
and the temporal evolution of image state electrons were performed using femtosecond time- 
resolved photoemission. A cluster of IPS with n=l is located near 0.75 eV below the vacuum 
level. These data are in agreement with results obtained in the course of studies of the 
emission properties of CNT |60l I66j . The values of the effective work function calculated 
from the Fowler-Nordheim plot are in the range 0.3 — 1.8 eV and the highest value of f3 FN 
extracted from the data is of the order of 1300. 
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3 The Electric Field of Relativistic Charges 



3.1 Single Charged Particle 

The electric field of a charged particle at rest is spherically symmetric. If the particle is 
moving with a uniform velocity v in an inertial frame S, its electric field is deformed: along 
the direction of motion the electric field becomes weaker by a factor 7 2 , while in the transverse 
direction the electric field is enhanced by a factor 7. Here, 7 denotes the particle Lorentz 
factor. 

Let (x,y,z) be the coordinates and r the position vector of the point P in S, the z-axis 
of the frame being taken along the direction of motion of the charged particle. Then the 
components of the electric field produced by a rapidly moving charge q are ^H], 

^ = r 2^ 2 1 K V ^2 13/2 ( X ' y> Z ~ Vt ^ ( 19 ) 

[x z + y + 7 [z — vty 

where the value of k depends on the system of units used, in our case k = 1 7 . We use 
the particle charge q rather than the elementary charge e to include particles with multiple 
charges such as ions for which q = Ze. At t = and z = rcosO, where 9 is the angle which 
the vector r makes with the z-axis, equation ()19|) can be rewritten in the form 



<?7 

K 2 



1 — (3 2 



3/2? - -> 

B ~ /3x F. (20) 



r 



.1 - (3 2 sin 2 9 

Thus, the magnitude of the electric field at 9 = 0° or 9 = 180°, Fn, is given by 



*ii = (21) 



while in the transverse direction, the magnitude of the electric field is 



F ± = (22) 



As seen from Eq. ()22|) . the increase of the strength of the electric field in the transverse 
direction corresponds efficiently to an increase of the particle electric charge, Q e ff = "jZe. 

In accelerators, charged particles move in front of conducting surfaces and one has to 
account for the image charge induced by the particle. In this section we discuss only the case 
when a positively charged particle is moving parallel to a grounded conducting plane. The 
particle and the image charge form a dipole. The complete form of the electric field vector 
of the relativistic dipole is presented in Appendix . 

Let the surface of the conductor coincide with the (yz)-plane and (0,y,z) be a point on 
the plane. The electric field at this point is normal to the surface and is directed into it. 
The components of the field from a positive point charge, located at the distance h above 
the plane, are (see Appendix) 

^ = — ^(-h,0,0). (23) 

Thus, at t=0 the magnitude of the electric field at the point directly beneath the positive 
charge is 

F, - ^. (24) 



Nevertheless, we keep the symbol k in equations below to control their dimensions. 
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We find that accounting for the image charge doubles the electric field strength on the surface. 

From equation (|24[) we conclude that only particles with ultra-high energies are able to 
produce significant field strengths at macroscopic distances. Let us substitute in (|24j) the 
numerical values of constants. That gives 

F x = 2.88 • 10~ 7 |j , (25) 

where h is the transverse distance in cm, F is the electric field in V/cm. For instance, the field 
strength of 1 V/cm at a distance 1 cm from the particle is generated by a positron of energy 
1.75 TeV or a proton of energy 3.25 • 10 3 TeV. Thus, only particles of ultra-high energies 
(e.g. cosmic rays) are capable to generate a large field strength at macroscopic distances. 



3.2 Bunch of Charged Particles 



Let us consider a bunch of TV" positively charged particles uniformly distributed within a 
circular cylinder of length L. Suppose that the bunch axis is along the coordinate z-axis and 
the bunch is moving along the z-axis with a relativistic velocity v = c(3. The radial electric 
self-field of such a rapidly moving bunch is described by |18j 



F±{r, z) 



rL l y / r 2 +J 



2^2 



1 + ^C 2 

1+ 8r 2 C l 



+ 



L 



s/r 2 +7 2 (i - zf 



3tf 



with 



Ci 



1 + 



2 2 

7 z 



-i 



C 2 



1 + 



l\L-zf 



-i 



(26) 



(27) 



The field described by (|26|). has a different behavior at distances far apart from the bunch 
and in the near region, r < L. At very large distances, r > 72 and r 3> j(L — z), equation 
(|26|) reduces to the Coulomb form (|22[). At the same time, in the near region and beyond 
the bunch tails, 72; j(L — z) 3> r and equation (|26|) simplifies to 

2qN 1 



L 



(28) 



which coincides with the external field of a continuous beam with the linear charge density 
A = qN/L. 

In an accelerator, a charged beam is influenced by its an environment (beam pipe, mag- 
nets, collimators, etc.), and a high-intensity bunch induces surface (image) charges or currents 
into this environment. This modifies the electric and magnetic fields around the bunch. For 
instance, if the bunch is moving in the midplane between infinitely wide conducting planes 
at x = ±/t (this simple geometry models a collimator gap), then the transverse component 
of the bunch electric field is described by |18j 

2qN it/2 



F±,tot{x) 



Lh 



(29) 



sin(f 6) ' 

where 5 = x/h. Thus, on the surface, S = 1, the field is enhanced by a factor tt/2 due to 
the presence of the conducting planes. If the bunch is displaced in the horizontal plane by x 
from the midplane the resulting field behaves as HF 



F± t tot(x,x) 



2qN 
J ~Lh 



cos(U) 



5-5 



(30) 



Bin(§tf)-sin(!*) 

Here 5 = x/h. Equations (|28 [) -(|30 [) demonstrates that although in the relativistic limit the 
electric field does not depend on 7, the field strength may be very high if iV is large and the 
product Lh is small. 
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Figure 2: View from the close distance to the LHC collimator jaw surface |79j . 

4 Field Emission in HEP Colliders 

In Section 2 the conditions were described under which FE starts to develop. Let us examine 
in this section the performance of the Large Hadron Collider (LHC) and the US Linear 
Collider (USLC) 67 , in an attempt to display their "bottlenecks". In the following, we 
shall not consider well established locations such as RF superconductive cavities, where 
FE routinely occurs and leads to the rapid growth of power dissipation with field. Our 
consideration here is limited to locations where bunches are moving close to conducting 
surfaces and create electric fields strong enough to activate observable FE effects. For the 
analysis we need to know the beam parameters, distances to the closest surfaces, their shapes, 
and materials their constituent. 

4.1 The LHC Collimator System 

Each of the two LHC rings will store 2808 bunches, each bunch populated with N = 1.15- 10 11 
protons at energies of up to 7 TeV (nominal design parameters). To handle the high inten- 
sity LHC beams and the associated high loss rates of protons, a strongly acting collimator 
system was designed and constructed |ZZ]-[HS]- The LHC collimators must restrict the me- 
chanical aperture in the machine and clean the primary halo so that quenches of magnets 
are avoided. To provide the required beam quality at the interaction points, primary and 
secondary collimators must be physically close to the high intensity beams. This impose 
strong limitations on the materials of the jaws and the length of the jaws. In the "final" 
solution for the collimation system [HS], |84| . the collimators in IR3 and IR7 (and IR6) are 
made of graphite with parallel jaw surfaces and small gaps, ~2 mm, at 7 TeV. Therefore, 
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the first candidate for a bottleneck is the very narrow collimator gap. 

Different scenarios for the beam parameters are foreseen in the LHC operation. These 
parameters are listed in Table 1 [BE]- For the nominal scenario, in accordance with (|29|) . the 
field strength at a distance of 1 mm from a bunch is equal 69 KV/cm. At this field level FEC 
is negligible if calculated with use of the classical FN equation l|4"|). However, the surface of a 
jaw made of graphite may be polished only mechanically and the resulting surface roughness 
is high (see Fig [2j) . The surface roughness and many other surface effects as discussed in 
Section 2 are characterized by the field enhancement parameter (3 FN , whose actual value 
needs to be determined by direct measurements. Thereby, the calculations of FEC presented 
below are performed for a range of (3 FN values. The FEC density, J FN ((3 FN F,T), calculated 
for the graphite jaws with use of (|29j) and (|12|) at T=300 K, h=0.1 cm and (/?=4.6 eV is shown 
in Fig. 1. For (3 FN below 300 and the nominal scenario (LHC-n), the electron emission is 
negligible. FEC for (3 FN in the range 300-500 begins to be notable and at higher (3 FN may 
became a very serious problem. In the advanced scenarios (LHC-0/1), when bunches are 
twice as short and their population is higher by the factor 1.48 (the field level on the surface 
increases by the factor 2.94), the emission current at /3 FJV =200 increases by 14 orders of 
magnitude ! At higher P FN , a very dense electron flow will cause electrical breakdown in 
vacuum, jaw surface heating and damage. The electron flow could also disturb the proton 
beam trajectory and give rise to loss of protons. The current density depicted in Fig. 1 
can be translated into the number emitted electrons, N e . A bunch with parameters from 
Table 1, expose the area 6A = 2L ■ a* during the time r = L/c will extract from both jaw 
surfaces N e = 25AtJ electrons. That gives N e = 7.93 • 10 7 J for LHC-n and N e = 1.36 • 10 7 J 
for LHC-0/1, respectively. Fig. 2 shows the evolution of N e with f3 FN - We note that at 
[i FN ~520 the number of extracted electrons N e (LHC-n) is comparable with the number of 
protons per bunch. 

Table 1 



LHC Beam Scenarios 85 , 86 , 87 


bunch parameters 


Nominal 


Phase 


Phase 1 


" Super-bunch" 


bunch population, N 
bunch radius, a* 
bunch length, L 
beam energy 


1.15xl0 u 
16.7 /im 
7.55 cm 
7 TeV 


1.7xlO n 
11.3 (im 
3.78 cm 

7.45 TeV 


1.7xlO n 

3.78 cm 
7.45 TeV 


5.6 xlO 14 

7500 cm 
7 TeV 



The FE current rises not only with F e f f , but the number of active emitters also increases 
with the field. A local increase of the emission current causes Joule heating of the emission 
tip which in turn leads to an increase of FEC. Fig. 3 illustrates the effect. It shows the 
temperature dependent factor in (|T2")) . whose rise is limited by the conditions (jTl)) and (fH)|) . 
For a given field strength and higher temperatures it is necessary to apply equations of the 
thermo-field-emission theory [5J E] • 

Joule heating by FEC is not the only mechanism responsible for an abrupt increase of 
FEC. Many studies point out that other mechanisms are also responsible for melting or 
damage of the emitter, such as the bombardment of the surface by ions. These ions are 
created by FEC as it passes through a cloud of gas evolving from the emission site |26( I27j . 
The plasma that develops can create an electric field on the order of 10 7 V/cm at its interface 
with the conducting wall and may result in an explosion of the emitter |28|-|31|.|20|. 

In Section 2 we noted that the image potential states on the surface of solids are charac- 
terized by low values of the work function, (p < 1 eV. If we hypothesize that the field emission 
due to IPS contributes in an important way to FEC, then the inner beam screen, consist- 
ing of copper tubes lining the LHC vacuum chamber, is another candidate for a bottleneck. 
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Indeed, the proton beam is moving at a distance of ~2 cm away from the beam screen wall 
with cryopumping slots |91| I92| , whose sharp edges could be very effective emission sites with 
high (3 FN - Fig. 4 shows the dependence of the FEC density on f3 FN calculated with 99=0. 92 
eV, /i=2cm and T=20 K, the temperature at which the beam screen will be maintained. As 
seen from the figure, for the LHC-n scenario IPS electrons will be seen only if (3 FN >500. 
However, in the LHC-0/1 scenario emission will be large already if (3 FN ~150 or higher. At 
(3 FN ki 600 the top of the surface potential barrier is located below the energy level of IPS 
with n=l. These electrons escape the surface freely and this is reflected as a change in the 
slope of the curve. In FigEJalso shown calculations for the super-bunch scenario (see Table 1) 
with the nominal work function for cooper, ip= 4.65 eV. 

The results presented in Fig El are to some extent speculative. The calculations were 
performed assuming that the density of IPS electrons, pumped by the soft component of the 
synchrotron radiation from the bulk to the surface is of the same order as the density of 
electrons at the Fermi level. Nevertheless they give rise to a cause for concern that FE from 
the beam screen may amplify electron cloud formation. 

4.2 The US Linear Collider 

In the similar way we may evaluate the level of the field emission in other colliders, for example 
proposed US superconducting linear collider (USLC) 67 which can became a prototype of 
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a future international linear collider (ILC). The USLC will bring into collision electrons 
and positrons at energy -^=500 GeV or, after an energy upgrade, 1000 GeV. The vacuum 
chamber in the positron ring will be coated throughout with a material with low secondary 
electron yield (e.g. conditioned titanium nitride, TiN) to prevent build-up of electron cloud. 
The beam pipe has a circular internal cross-section, with radius h=21 mm. The vacuum 
chambers are all constructed from an aluminum alloy. The bunch length is extraordinary 
short. The bunch compressors must reduce the ~5 mm rms length of the bunches extracted 
from the damping rings to 300 |im bunch length required for the main linac and final focus 
systems. For optimum collider performance in the cold option, bunches with a charge of 
N = 2.0 • 10 10 and a length of L = 300/mi will be injected for acceleration in the main linac. 

For a beam well centered in the circular beam pipe, the image charge effect is vanished 
[55] . For this reason, instead of $Z§§ we use in calculations equation (|2*5)) . As a result, the 
electric field strength on an ideal smooth surface is characterized by the parameter N/hL. 
Direct calculations shows that 



denoting that with the same (3 FN and ip as for the LHC, the field emission from the USLC 
beam pipe surface will be stronger than in the LHC collimators. Fig. shows the FE current 
density predicted for a variety of <p values. The range of (p employed in the calculations reflects 
uncertainties in the work function of TiN reported by different authors ^6Q|-[Z£]. There is an 
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200 400 600 800 1000 1200 

T,K 

Figure 5: The temperature increment of the FE current (|12|) at different f3 FN in the temperature 
range as restricted by the conditions (|14 |l -(|15 j> . To denotes the temperature limit as determined 
from JI3J). 



observation [TUj, that changing the nitrogen concentration can modify the work function of 
TiN. Furthermore, a thickness dependency of the work function was observed for TiN/Al 71 . 
For a TiN layer with thickness less than 20 angstroms, the TiN/Al gate work function is the 
same as reported Al work- function (about 4.08 eV). For TiN layers of thickness greater than 
100 angstroms, the TiN/Al gate work function is approximately the same as reported TiN 
work-function (4.5 V). For TiN layers with thickness from 20 angstroms to 100 angstroms, 
the work function can change as the TiN thickness is changed. The TiN crystal orientation is 
important because the metal work function is strongly dependent on this quantity. According 
to [Z2]) (100) orientation TiN has <p ~4.6 eV and 4.4 eV in the (111) orientation. The work 
function of the TiN film was found to be sensitive to the film morphology, stipulated by a 
sputtering technology and may be as low as 3.72 eV [73] . 

In this way, the uncertainty in the value of (p gives rise to an uncertainty in the current 
density of 3-5 orders of magnitude. 

We finally point out that parameters of the positron and proton bunches and the beam 
pipe geometry at the HERA collider are such that even at (5 FN = 10 3 the emission current 
density is well below 10~ 17 A/cm 2 . 

5 Electron Dynamics in the Bunch Fields 

In this section we discuss in a semiquantitative way the fate of electrons after their emission 
in the LHC collimator gap. Electrons leaves the surface with different energies. The energy 
spectra of electrons just ejected by the surface is presented in Fig El For a given imposed 
field, the width of the energy distribution grows with an increase of the surface temperature 
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Figure 6: Densities of the electron current ejected from the copper beam screen at T=20 K: 
(LHC-n, LHC-0/1) due to presence of the surface image-potential states (IPS), <^(IPS)= 0.92 eV, 
LHC-sb) the super-bunch scenario , <p(Cu)= 4.65 eV. 



and more energetic electrons easier escapes the emitter. 

Let us consider a single electron motion in the self-fields of a bunch. The analysis is based 
on the classical Lorentz equation [SSI EH ISO] 

|= eF + ec^xB], (32) 

To be specific, assume that the electron is emitted from the plane surface with the velocity 
vector, = v^ e ^/c, normal to the surface. The surface is placed parallel to the (0, y, z) 
plane at a distance x = +h. The bunch is moving along the line (0,0, z) in the positive z 
direction, as shown in Fig. |SJ In the simple geometry chosen, we consider, only as an illus- 
tration, the evolution of an electron trajectory in the (x, 0, z) plane; in this way components 
of the self-field vectors can be set F = (F x ,0, 0) and B = (0,B y ,0), where F x is of the form 
I® and dl 

' -/?j 6 >. (33) 



_ *2qN 
y ~ c Lh 



B, , s tt/2 

-£ + (l-BS 



Here fti^ is the z-component of the bunch velocity, 5 = x/h, B. = ^/2ttL/2t:R is the 
Laslett bunching factor and R the average LHC machine radius. In the LHC-n scenario, 
B, = 7.05 x 10~ 6 HHIIM] and terms proportional to B. can be omitted in (|3*3*|) . 

Surface irregularities of a very small size possess tips that cause a local field enhancement 
and increased emission. The field of such a tip decays with distance approximately as [TU] 

,p + h\l-5\. 



M(5) = 1 + P FN ( - ■ P jr ) , (34) 
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Figure 7: Densities of the electron current emitted from the USLC main linac beam pipe at T=20 K 
and various values of the work function. 



where S = x/h is a normalized distance from the emission tip, p the tip radius. The self- fields 
F and B enter equation (|3*2*|) multiplied by the factor • 

For numerical calculations, it is convenient rewrite (|32j) by use the dimensionless time 
r = t/rQ with To = L/c. For a full generality, we retain a relativistic treatment of equation 
()32|) which in components becomes 



d(/#7 (e) ) 
dr 

d(M e) 7 (e) ) 
dr 



where the function -F(<5) outside the bunch takes the form 

M(S) 



F{5) = A ■ 



sm(UY 



a < \x < h, 



and inside the bunch is H 



with the dimensionless constant 



2A ( h 2 \ . 



\x\ < a 



kit eqN 

A = pr-r . 

mc La 



(35) 



(36) 



(37) 
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A direct estimation gives that A = 0(1) for the LHC-n parameters. Fig El schematically 
shows electron trajectories and forces acting on the electron as a consequence of equation 
©■ 

Let us analyze the system (|H5[). For a relativistic bunch one may put fig = +1 and 
therefore below we omit superscripts in the electron velocity components. A partial solution 
of (j.35j) can be found for nonrelativistic electrons, 7 ~ 1, by taking a ration of the two 
equations. That gives 

^ = (0,-1)^. 
ar ar 

After integration of the last equation and some algebra one finds a solution in the form 



0z = 1- V 1 + Pl~Ph 

where the integration constant (5q is fixed by initial conditions (5 Z =0, /3 X = 0o a t t=0. This 
equation connects the p z and p x components of the electron momentum in the following way 

» = -4s£ (38) 

Equation (|88jl implies that during the electron acceleration p z remains negative and much 
smaller than p x . Thus, the electron displacement in the z direction is expected to be small 
too. 

The magnetic field of the bunch gives rise to the p z component and the electron turn in the 
(x,0,z) plane. For an accelerated charge the radiation emitted at any instant is approximately 
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Figure 9: Illustration to Eq. (|32|) . Electrons ejected from the surface with the velocity v = (v x ,0,0) 
are accelerated in the radial electric field and rotated by the magnetic field B, created by a cylindrical 
proton bunch. The Lorentz force acting on the electron is denoted by F L and the beam current by Jt>- 



the same as that emitted by a particle moving instantaneously along the arc of a circular 
path whose radius of curvature is Ri. The frequency spectrum of radiation emitted by a 
charge in instantaneously circular motion is characterized by the critical frequency 

^ c = ^7 3 , (39) 

beyond which there is negligible radiation at any angle |88| I9()j. The corresponding energy 
of the photon is E = hu> c . 

For a nonrelativistic charge whose motion is described by the system l[35JI with 7^) = 1, 
the instantaneous radius Ri can be found in the following way. Let the point (xq, Zq) be the 
center of the circle of radius Ri 

(x - x Q f + (z- zof = R\ (40) 

whose infinitesimal arc with the point (x, z) on it coincide at the moment r with the electron 
trajectory. To find Ri we differentiate (|%U|) twice by r and find relations between x — xo, 
z — zq and f3 x , f3 z , (3 X and (3 Z . Using (|35]> for (3 X and (3 Z , we obtain 

Rl{T) = WM^WTpTY (41) 

Now we are in a position to present numerical results. The input parameters for a 
numerical solution of the system IjMftjl are listed in Table 2. These are the distance, h, 
from the bunch center to the collimator plane, the tip radius p in equation (|34j). the field 
enhancement factor at the surface, (3 FN , the electron velocity (3, or its kinetic energy E kin . 
Parameters of the proton beam are correspond to the LHC-n scenario, see Table 1. 
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Table 2 



Input paramerets for the system (j35j) 


h, mm 


P, m 




. , eV 






1 


1.0xl0~ s 


300 


4.6 


-4.2xl0- 3 


0.0 



Table 3 presents parameters of the electron trajectory, whose numerical values were cal- 
culated at the emission point (5=1), at the bunch center (5=0) and at the opposite collimator 
plate (5 = —1). The trajectory is split into three parts. At the first stage, the electron is 
accelerating by the field (|36|) and moves toward the bunch (FigEJJ). The local field enhance- 
ment (|34p causes an extremely high acceleration rate and a rapid increase of the electron 
momentum. At the start time the instant curvature of the trajectory Ri is very small and 
the energy loss E through radiation can run up to 5% of the initial electron energy. At later 
times with an increase in the distance from the surface, the energy loss by radiation is very 
low, E 1 ~ 10~ 6 - 10- 4 eV. 

At the second stage, the electron is crossing the bunch. The internal field of the bunch is 
described by (|37|) . At the bunch center the electron momentum reaches a maximum value of 
151 KeV/c at r =0.129 after emission. An average distance between protons in the LHC-n 
bunch is of the order of 100 nm and the electron either passes through the bunch or will 

Table 3 



Selected paramerets of the electron trajectory 



5 


T 


P 


E, - , eV 

kin * 


p, KeV/c 


Ri, m 


iv c, Hz 


E„ eV 


■s/sepi GeV 


1 





4.2xl0- 3 


4.6 


2.17 


l.lxlO" 6 


4.3xl0 14 


0.283 







0.129 


0.293 


2.3xl0 4 


151.0 


7.7xl0~ 2 


6.7xl0 9 


4.4xl0~ 6 


2.94 


-1 


0.254 


3.3xl0~ 2 


271.1 


16.6 


4.8xl0~ 3 


9.5xl0 10 


6.2xl0~ 5 





be scattered by a proton. In the later case, the electron and the proton undergo an elastic or 
inelastic collision at the center of mass energy of 2.94 GeV. In the rest frame of the proton, 
the electron carries the momentum of 4.13 GeV/c. Such an amount of energy is enough for 
hadron production in an inelastic ep collision. 

At the last stage, the electron is decelerated by the bunch field. If the electron passes the 
bunch without scattering, then at the moment when r =0.254 it impacts the opposite plate 
of the collimator at 5 =-1 with a kinetic energy of 271 eV. For a fully symmetric configuration 
and neglecting the electron energy loss in the course of acceleration and deceleration, the 
electron has to arrive at the opposite plate with the same kinetic energy as at the start point. 
The energy difference arise due to an additional acceleration of the electron in the enhanced 
field at the tip apex and the assumption that at the end point the surface is flat. 

Above, we considered the emission of a single electron, when at r=0 the emission point and 
the bunch head had the same z-coordinate. However, the field emission occurs throughout 
the entire region L, covered by the field of the bunch. Each electron reaches the bunch 
center after time At ~ 0.129. If now account a bunch movement, then electrons emitted at 
1 — 2 At < t < 1, will experience only initial acceleration or deceleration in the wake field of 
the bunch. Thus, the electron energy spectrum will be very broad, from 4.6 eV to 23 KeV. 
As a result of emission from both collimator plates, the proton bunch occurs in the electron 
streams being intersected. For high f3 FN it is necessary take into account the space charge 
of electron clouds, whose electric field reduces the the strength of the bunch field and leads 
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to a decrease of the electron emission. Nevertheless, the leading part of the bunch of length 
At will always remains out of the electron cloud and triggers the field emission. 

The interaction between the electron and the surface results in the ejection of secondary 
electrons from the material. The secondary electrons consist of true secondaries and those 
elastically reflected. The number of secondary electrons is given by the secondary emission 
coefficient (SEC) that depends on the surface characteristics and on the impact energy of the 
primary. In turn, the secondaries are accelerated and, on impact, produce further generation 
of electrons. The electron-cloud build-up is sensitive to the intensity, spacing, and length 
of the proton bunches. The electron flow increases exponentially if the number of emitted 
electrons exceeds the number of impacting electrons, and if their trajectories satisfy some 
specific conditions. For most materials the SEC exceeds unity for impact energies in a 
range from a few tens of eV to a few thousand of eV. This type of electron multiplication, 
so-called multipacting, and build-up a cloud of electrons has been studied recently very 
intensively |91j- jl()4] . In the LHC arcs, as believed, the dominant source of electrons will be 
photo-electrons from synchrotron radiation and beam-induced multipacting to be the leading 
source of sustained electron-production. 

The region near the scrapers and collimators is susceptible to a high beam-loss, and is 
potentially another location of high electron concentration. Protons incident on the colli- 
mator surfaces produce secondary electrons. Depending on the energy of the beam and the 
incident angle, the secondary electron-to-proton yield can greatly exceed unity when the 
incident beam is nearly parallel to the surface. 

In conclusion let us note that any heterogeneous structures such as a grooved metal 
surface |1()5[ IT06_ introduced with intention to combat with the multiplication of electrons and 
located near an intensive positively charged beam, may in real conditions actually provoke 
field emission. 

6 Praemonitus Praemunitus 

The analysis performed in the previous sections makes it possible to draw the conclusion 
that in the addition to the known electron sources, electron field emission intensified by 
multipacting can make a dominant contribution to the build-up of an electron cloud in 
the LHC collimator system and may become a serious problem. The analysis of the ILC 
prototype reveals that a noticeable field emission will accompany positron bunches on their 
entire path during acceleration. 

From the examples considered we learn that the level of field emission is controlled by 
two essential parameters. The first of these, (3 FN , is wholly determined by the state of the 
emitting surface and by the processes proceeding on it. Therefore, control of f3 FN is possible 
only within certain limits. The effect of surface aging will increase (3 FN with time. The actual 
value of (3 FN for the graphite jaw can be obtained only by direct measurement. 

For P FN below 300 and the nominal LHC-n scenario, the electron emission is negligible. 
If the j3 FN value is in the range of 300-500, FEC begins to be notable. At higher (3 FN , FEC 
may became a very serious problem. In the advanced scenarios (LHC-0/1), when bunches 
are twice as short and their population is higher by a factor 1.48, the emission current at 
/3 FJV =200 increases by 14 orders of magnitude as compared with LHC-n. At higher f3 FN , a 
very dense electron flow will cause electrical breakdown in vacuum, jaw surface heating and 
damage. The electron flow could also disturb the proton beam trajectory and give rise to a 
loss of protons. 

The value of the second parameter, N/hL, is assigned by the beam parameters and the 
distance to the conducting surface. It is possible and necessary to choose the components of 
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/? FJV N/hL in the way to minimize the electron field emission. 
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Appendix: Electric Field of a Relativistic Dipole 



Take two point charges, +q and -q, separated by a distance 2h. The charges are at rest in 
frame S and the x-axis passes through the charges. We denote the coordinates of a point in 
the frame S by (x,y,z). Then the potential from the two charges (a dipole) is given by 



y/[x - h) 2 + y 2 + z 2 y/[x + h) 2 + y 2 +z 2 
We obtain the electric field of a dipole at rest by the usual rule: 



}■ 



or in components 

x — h x + h 



F x = Kq 
where 



D 3 _ 



F y = Kq 



1 1 



y, F z = Kq 



r l l 



D± = y/(x ± h) 2 + y 2 + z 2 . 



(A-l) 



(A-2) 



(A-3) 



(A-4) 



Now, suppose that the dipole is moving along the z-axis of the frame S with velocity v. Then 
the space coordinates in S and S are related by the Lorentz equations 

x = x , y = y, z = j(z-vt), (A-5) 

The transformation laws for the components of the electric field can be written as 

F x = jF x , F y = jF v , F z = F z . (A-6) 

Substituting expressions of F from (|A-3|) into ()A-6j) and account (|A-5[) . we arrive at 

'x — h x + h~ 



Kq>y 
K Q1 



i i 

- 1 1 

.W ~ dJ 



Dl 

y 



where 

On the plane x = 



Fr 



(z - vt), 

D± = y / {x±h) 2 +y 2 +j 2 {z -vt) 2 . 

2Kqjh 
[h 2 + y 2 + 7 2 (z - vt) 2 ] 



2 13/2' 



F — F 



0. 



(A-7) 
(A-8) 
(A-9) 

(A-10) 
(A-ll) 
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